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Summary 

Nonlinear mathematical programming methods 
are used to design a radiantly cooled and heat-pipe- 
cooled panel for a Mach 6.7 transport. The cooled 
portion of the panel is a hybrid heat-pipe/actively 
cooled design which uses heat pipes to transport 
the absorbed heat to the ends of the panel where 
it is removed by active cooling. The panels are op- 
timized for minimum mass and to satisfy a set of 
heat-pipe, structural, geometric, and minimum-gage 
constraints. Two panel concepts are investigated: 
cylindrical heat pipes embedded in a honeycomb 
core and an integrated design which uses a web-core 
heat-pipe sandwich concept. The latter concept was 
lighter and resulted in a design which was less than 
10-percent heavier than an all actively cooled con- 
cept. The heat-pipe concept, however, is redundant 
and can sustain a single-point failure, whereas the 
actively cooled concept cannot. An additional study 
was conducted to determine the optimum number 
of coolant manifolds per panel for a minimum-mass 
design. 

Introduction 

The use of nonlinear mathematical programming 
methods (formal optimization) for the synthesis of 
aeronautical and/or structural systems is fairly well 
established (refs. 1 and 2), but the use of such meth- 
ods to design fluid-loop and/or heat-pipe systems 
has been limited (e.g., ref. 3). This paper is in- 
tended to demonstrate an efficient means for design- 
ing optimum, multifunction, heat-pipe systems us- 
ing formal optimization methods, with the hope that 
such methods can replace cumbersome parametric 
methods (refs. 4 and 5). 

The application chosen for study is the use of heat 
pipes to cool structural panels of a hypersonic cruise 
vehicle. This particular heat-pipe application was 
chosen for several reasons: it incorporates several 
disciplines of engineering analysis (thermodynamics, 
fluid dynamics, and structures), it is a viable applica- 
tion of heat pipes for future hypersonic vehicles, and 
it demonstrates a method for designing heat pipes to 
perform several functions efficiently (e.g., transport 
heat and support structural loads). 

References 6 and 7 describe a combined radi- 
antly and actively cooled panel (RACP) applicable 
to a Mach 6.7 hypersonic transport. The system 
uses heat shields and insulation to substantially re- 
duce the incident heat load to the cooled structure 
and results in an overall structural design which is 
safer, more tolerant to off-design conditions, slightly 
lighter, and easier to fabricate than a design which 
must absorb all the incident heating. (See ref. 8.) 


However, the cooled structure remains susceptible to 
damage or leakage problems, which requires a redun- 
dant system, and hence additional weight and com- 
plexity, to avoid jeopardizing the flight of the ve- 
hicle. This paper is an investigation of the design 
of a heat-pipe-cooled panel which could replace the 
cooled-panel portion of the RACP described in refer- 
ences 6 and 7 (i.e., that portion of the panel between 
coolant inlet and outlet manifolds). An important 
advantage of using heat pipes for such an application 
is the increased redundancy in the design (e.g., if 
one heat pipe fails, the entire panel/vehicle will not 
be lost) with less weight penalty and system com- 
plexity. Heat pipes do not contain large quantities of 
working fluid, and if one should rupture, fluid leakage 
would cause minimal damage. Also, heat pipes natu- 
rally tend to eliminate any thermal gradients/stresses 
which might occur because of high localized heating. 

Designs are obtained for several heat-pipe-cooled 
panels (HPCP) using formal optimization methods. 
In each design, heat pipes are used to transport the 
heat load experienced by the cooled panel to coolant 
manifolds located along the panel ends. The coolant 
manifolds for the heat-pipe-cooled panels contain 
both supply and return lines and circulate fluid and 
heat back to heat exchangers located on the vehi- 
cle. The effects of this redesigned manifold on to- 
tal system mass and structural performance were 
not considered, and no attempt was made to in- 
clude the coolant manifold design or coolant circu- 
lating system in the optimization. The panels are 
designed for minimum mass and to satisfy a set of 
heat-pipe, structural, geometric, and minimum-gage 
constraints. The first concept consists of cylindri- 
cal heat pipes embedded in a honeycomb sandwich 
structure. This design represents a retrofit, where 
the cylindrical heat pipes would replace the D-tubes 
of the previous RACP design. (See ref. 7.) It is 
not a truly integrated design (thermally and struc- 
turally), because the heat pipes, although capable 
of sustaining structural load, remain self-contained 
units which are bonded to the structure as opposed 
to being an integral part of the structure. Also, for 
this design the heat pipes are sized to satisfy heat- 
transport and fatigue requirements, and the sand- 
wich panel is then resized to accommodate the struc- 
tural loading. The fatigue requirement is satisfied 
by requiring the heat-pipe tube wall to be thicker 
than some minimum value. The second concept is 
a web-core heat-pipe sandwich panel which is opti- 
mized for minimum mass while subject to heat-pipe, 
strength, buckling, fatigue, and minimum-gage con- 
straints. This concept represents an integrated heat- 
pipe/structural design where the heat pipe is an in- 
tegral part of the structure and is designed for two 


functions simultaneously — to transport heat and to t 
support structural loads. Both panels were designed 
to withstand the same set of thermal/ structural load- M 

ing conditions and constraints as the RACP of ref- 
erence 7. Further, the effect of additional, equally M x?max 
spaced coolant manifolds on optimized heat-pipe- 
panel mass was studied. The design methodology 
used in the present study to design the heat-pipe m 

panel is compared with that presented in reference 7 iV 

for designing the RACP. N 

The optimizer used in the present study is ,n 

CONMIN (ref. 9), which uses the method of feasible- 
directions algorithm for constrained function mini- p 
mization. The optimizer is coupled to analysis and 
plotting modules via FRANNOP (a computer pro 


gram described in ref. 10). 

J X 

Symbols 


Px,ob 

A 

cross-sectional area, in 2 


a 

crack length, in. 

Pv 

do 

sonic velocity, in/sec 

Q 

b 

panel width, in. 

Qain . Qtin 

D 

groove depth, in. 


D x 

panel bending stiffness, in-lb 

QawiQasi 

Qae 

da/dN 

crack growth rate, in/cycle 

E 

modulus of elasticity, lb/in 2 

Qin 

F 

objective function, lbm/in 2 

Qtw 

FS 

factor of safety 

/maxi /min 

maximum and minimum cyclic 

Qx, max 


stresses, ksi 

9c 

gravitational constant, 


386 lbm-in/lbf-sec 2 

q 

9i 

ith constraint 


h 

heat-transfer coefficient, 
Btu/ft 2 -hr-°F 

Qabs 

he 

panel-core height, in. 

R 

lx 

moment of inertia about z-axis 

o 

per unit width, in° 

ri,r2,r 3 

Kc 

stress-intensity factor, ksi-\/im 

r p 

K w 

wick permeability, in 

r v 

AK 

stress-intensity range, ksi-\/im 

q 

L 

panel length, in. 

o 

T 

£ e ff 

effective length of heat pipe, 
in. 

I 

t 


simply supported plate length, 
in. 

mass, Ibm 

maximum panel bending 
moment per unit width, in * n *^ 

number of half waves 

number of load cycles 

number of additional coolant 
manifolds 

spacing between heat pipes or 
coolant tubes, in. 

applied in-plane load in x- 
direction, lbf/in. 

overall panel buckling load per 
unit width, lbf/in. 

rt 

internal pressure, lbf/in 
flaw shape parameter 

axial and transverse heat input 
loads, Btu/sec 

axial heat-pipe wicking, sonic, 
and entrainment limits, 

Btu/sec 

input heat load to panel, 
Btu/sec 

transverse heat-pipe wicking 
limit, Btu/sec 

maximum panel shear load 
per unit width in x-direction, 
lbf/in. 

uniform transverse panel load 
per unit area, lbf/in 2 

absorbed heat input to cooled 
panel, Btu/ft 2 -sec 

stress ratio 

heat-pipe radii shown in 
figure 6, in. 

pore radius of wick, in. 

hydraulic radius of heat pipe, 
in. 

spacing of panel-core webs, in. 
temperature, °F 
plate thickness, in. 
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tf i tc 

face-sheet and core thick- 
nesses, in. 

to 

outer facing thickness, in. 

tsi tyj 

screen and wick thicknesses, 
in. 

V 

static moment of area, in 3 

w 

groove width, in. 

w 

simply supported plate width, 
in. 

Xi 

ith design variable 

y 

distance from neutral axis, in. 


coefficients 

7 

ratio of specific heats 

A 

latent heat of vaporization, 
Btu/lbm 

P 

viscosity, lbm/sec-in. 

V 

Poisson’s ratio 

P 

mass density, lbm/in 3 

<?e 

liquid surface tension, lbf/in. 

^FATi a x,p, & Xi lb, 

stresses defined in structural 

max, 0y,p, CTfc c , 
^x,/ 1 Tzx,ma.x 

analysis, lbf/in 2 

^ult > r ult 

ultimate normal and shear 
stresses, lb/in 2 

* 0 

Subscripts: 

buckling parameter 

c 

core (web) 

f 

face sheet 

£ 

liquid 

lb 

local buckling of face sheet or 
webs 

ob 

overall panel buckling 

ref 

reference 

s 

solid 

V 

vapor 

w 

wick 


Design Problem Description 

A radiantly and actively cooled panel (RACP) 
was designed for a location 10 ft aft of the nose on 


the lower fuselage centerline of a Mach 6.7 hydrogen- 
fueled hypersonic transport described in reference 7. 
The design requirements and criteria for the RACP 
of reference 7 are discussed subsequently with a 
description of design methodologies for an RACP 
design and an alternative version of it, a radiantly 
and heat-pipe-cooled panel (RHPCP), which is a 
hybrid design using active cooling and heat pipes to 
cool the structural panel. 

Design Loads and Criteria 

Reference 7 describes the design requirements, 
design criteria, and design process for a full-scale 
2- by 20-ft RACP, shown in figure 1. The high-heat- 
Hux panel of reference 8 was designed for an incident 
heat flux of 12 Btu/ft 2 -sec to a 300°F surface tem- 
perature. A heat-transfer coefficient of 16 Btu/ft 2 - 
hr- F and an adiabatic wall temperature of 3000° F 
were used in reference 7 to define the heating envi- 
ronment. These quantities were used to determine 
heat input to the surface of the RACP and heat ab- 
sorbed q aba by the actively cooled-panel (ACP) por- 
tion of the RACP. It was determined in reference 7 
that a ^minimum- mass design results for q abs = 0.8 
Btu/ft -sec. A uniform surface pressure of '±1.0 
P®' a ” d r/ a c y clic . uniform, in-plane limit load of 
1200 lbf/in. along the 2-ft edge comprised the struc- 
tural loads on the panel. Design life for the RACP 
called for a 10000-hr exposure to operational tem- 
peratures, 300°F for the aluminum ACP, and 20000 
cycles (5000 cycles with a scatter factor of four) of de- 
sign limit loads. Stresses in the aluminum sandwich 
panel were restricted to values below which an ini- 
tial flaw size of 0.005 in. at the edge of a fastener 
hole would not grow to a critical length and sur- 
face flaws of 0.005 in. would not grow through the 
thickness of the coolant tubes or manifolds during 
the 20000 limit-load cycles. A factor of safety (FS) 
o 1.5 was applied to in-plane and transverse loads to 
obtain ultimate loads, and an FS of 4.0 was used to 
obtain ultimate burst-pressure loads. The radiantly 
and actively cooled panel was designed to sustain de- 
sign limit loads without yielding and design ultimate 
loads without structural failure. For the heat-pipe- 
cooled-panel (HPCP) portion of the RHPCP, an FS 
of 1.5 was also applied to thermal loads to ensure 
thermal redundancy if a heat pipe failed. 

Design Methodologies 

Two different design methodologies were used to 
design a cooled panel: parametric and trade studies 
were used to design the RACP of reference 7, and 
nonlinear mathematical programming (formal opti- 
mization) techniques were used to design the hy- 
brid actively cooled and heat-pipe-cooled panel of 
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the present study. Although the scope of the HPCP 
design study is limited, the study demonstrates the 
methodology as a viable alternative to parametric 
techniques, which are limited in effectively selecting 
an optimum design as the number of design variables 
increases. 

The parametric and trade studies used to design 
the RACP of reference 7 are illustrated in figure 2. 
Figure 2(a) indicates the overall process, and fig- 
ure 2(b) shows the numerous parametric and trade 
studies required to arrive at a final design. Because 
two mechanisms were used to accommodate the in- 
cident heating, i.e., rejection by radiation and ab- 
sorption by convection, it was necessary to establish 
the performance of each subsystem (heat shield and 
actively cooled panel). By using simplified analyses 
and the series of trade studies shown in figure 2(b), 
it was determined that a minimum-mass system re- 
sults when the cooled panel absorbs 0.8 Btu/ t '® ec - 
Once the major parameters and absorbed heat flux 
were determined, more extensive analyses and trade 
studies were performed to substantiate the design, 
to size the cooling system, and to arrive at the 
minimum-mass cooled structural panel which could 
absorb 0.8 Btu/ft 2 -sec. 

The current investigation replaces part of the 
trade studies with nonlinear mathematical program- 
ming techniques to design a heat-pipe-cooled panel 
to replace the portion of the cooled panel between 
the coolant manifolds. Although this represents only 
a portion of the total RHPCP system, it is suffi- 
cient to demonstrate the methodology. For this prob- 
lem, therefore, results of references 7 and 8 were 
used to determine the design and operating condi- 
tions. The absorbed heat load and maximum operat- 
ing temperature were assumed to be 0.8 Btu/ft -sec 
and 300°F, respectively. The optimization procedure 
is illustrated in figure 3. The procedure is discussed 
in detail in the section entitled “Statement of the 
Optimization Problem.” 


Actively Cooled Panel (ACP) Description 

The ACP, shown in figure 4, consists of an alu- 
minum honeycomb-sandwich structural panel with 
coolant tubes (D-tubes) next to the outer skin. The 
honeycomb structural panel is cooled by a 60-percent 
mass solution of ethylene glycol in water pumped 
through the D-shaped cooling tubes at a flow rate o 
3.4 gal/min. The panel is 2 by 20 ft and is supported 
by channel-section frames every 2 ft (representative 
of typical transport construction support) and J™ 3 
inlet and outlet coolant manifolds at the ends of the 
panel. 


Heat-Pipe-Cooled Panel (HPCP) Description 


The HPCP proposed in the present study differs 
from the ACP of reference 7 in two areas. Instead 
of an actively cooled structural panel (fig. 4), a pas- 
sive, self-contained, heat-pipe-cooled panel (HPC ) 
is used to transport heat to coolant manifolds on ei- 
ther end of the panel. Also, the manifolds of the 
HPCP include both supply and return coolant lines 
and circulate fluid and heat back to heat exchangers 
located on the vehicle. 

Heat pipes, shown schematically in figure 5, 
are self-contained devices which can transport heat 
nearly isothermally, at high rates, and over long dis- 
tances (in this case 10 ft) without the need for ex- 
ternal pumping. They operate by absorbing heat by 
the evaporation of a working fluid in a net heat inpu 
region, which in this case is the surface area of the 
panel. The heated vapor thus formed then flows to 
cooler condenser sections (actively cooled ends of the 
panel where coolant manifolds are located), where 
heat is rejected by the condensation of the vapor- 
The cycle is completed with the return flow of liquid 
condensate back to the heated region by the capi 
lary action of the wick. Several fluid-dynamic limits 
which can prevent normal heat-pipe operation are 
described in appendix A. 


The coolant manifolds on the HPCP would cir- 
culate the coolant along the ends of the panel, ex- 
tracting heat from the heat pipes, and back to heat 
exchangers located in the vehicle interior instead of 
through the panel as is done on the ACP. It is as- 
sumed that there is no mass penalty associated with 
the redesigned coolant manifolds. The heat-pipe 
panel offers more redundancy than the ACP design, 
because adjacent heat pipes are sized to accommo- 
date the additional heat load if a heat pipe fai s. 
Hence, loss of one heat pipe would not cause leakage 
of large amounts of coolant, catastrophic overheat- 
ing and subsequent total panel failure and a mission 
abort or vehicle loss. Two heat-pipe-panel designs 
are investigated: one design simply replaces the D- 
tube coolant passages of reference 7 with cylindri- 
cal heat pipes (fig. 6), and the other is an integral 
heat-pipe/structural configuration called a web-core 
heat-pipe sandwich panel (fig. 7). Two variations 
of wick configurations were investigated for each de- 
sign. Also, effects of adding more coolant manifolds, 
equally spaced along the panel, on the mass of the 




Cylindrical heat pipes embedded in a honeycomb 
sandwich. Cylindrical heat pipes were chosen instead 
of D-shaped heat pipes because the net heat load to 
the panel is small enough that temperature gradients 
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between adjacent heat pipes are of little concern and 
because the line contact of the cylindrical tube is of 
no consequence. Also, for the heat-pipe design, the 
cross-sectional area must be sufficient to accommo- 
date the axial heat transport, and the circular cross 
section offers more area than the D-shaped tubes. 
The honeycomb-sandwich dimensions were similar to 
those of the ACP, and the wall thickness of the heat 
pipes was maintained at 0.02 in. to satisfy fatigue 
constraints identified in reference 7. The heat pipes 
use a screen-covered, axial-grooved wick; the screen 
covering the grooves is a 165 x 1400 mesh screen with 
a pore radius r p of 0.0009055 in. (ref. 11). An alter- 
nate wick which consisted of screen-covered corru- 
gated channels where the channels and screen cov- 
ering use the same 165 x 1400 mesh screen was also 
investigated. The latter configuration (referred to 
as corrugated sandwich screen wick) increases the 
cross-sectional liquid area but decreases the amount 
of load-carrying material. The material is aluminum, 
and the working fluid is toluene. 

Heat-pipe web-core sandwich panel. The heat- 
pipe web-core sandwich panel design integrates the 
heat pipes with the structure to produce an effi- 
cient thermal/structural design. The web-core panel 
(fis- 7) consists of two face sheets of equal thickness 
tf and channel-section webs of thickness t c . Two 
wick concepts were investigated: screen-covered ax- 
ial grooves (165 X 1400 mesh screen, fig. 7(a)) and 
closed channels formed by using a sandwich wick 
of corrugated screen bonded between a layer of flat 
screen (165 x 1400 mesh) and the inner surface of 
the face sheets (fig. 7(b)). Both concepts assume 
a composite screen wick along the webs to trans- 
port fluid to the heated face of the panel. The com- 
posite screen wick (fig. 7) is composed of one layer 
of high-density screen (165 x 1400 mesh) along the 
liquid/vapor interface against several layers of 
coarser screen (325 x 325 mesh) along the web walls. 
The material is aluminum, and the working fluid is 
toluene. 


Statement of the Optimization Problem 

The design of heat-pipe panels is posed as a non- 
linear constrained mathematical optimization prob- 
lem, where the objective function to be minimized 
is the panel mass per unit planform area. The con- 
straints of the problem are the thermal, structural, 
geometric, and minimum gage requirements which 
the panel must satisfy. The design variables in each 
problem are selected cross-sectional dimensions of 
the panel. 


The mathematical statement of the optimization 
problem is as follows: find the set of design variables 
X such that 

F{X) minimum 
subject to 
9i(X)< 0 

where F is the objective function, in this case panel 
mass per unit planform area; where the design 
variables A are various cross-sectional dimensions 
of the panel; and where the inequality constraints 
gi(X) < 0 are various thermal/structural require- 
ments or limits. The optimization algorithm used 
in the present study is CONMIN (ref. 9), which is a 
general-purpose optimization program that performs 
constrained minimization using the method of feasi- 
ble directions search algorithm. Figure 3 is a flow 
chart of the optimization procedure. The procedure 
begins with the selection of an initial set of design 
variables. An analysis module performs the neces- 
sary thermal/structural analyses and the optimizer 
then calculates the objective function, constraints 
and the gradients of the objective function and con- 
straints with respect to design variables. Gradients 
required by the search algorithm are calculated by 
finite differences. Using the method of feasible direc- 
tions search algorithm, the optimizer selects a search 
direction and minimizes the objective function in this 
direction. The design variables are modified to reflect 
this design choice and the procedure is repeated until 
one of several termination criteria is satisfied The 
termination criteria are: (1) the percentage change 
in the objective function is less than some specified 
value, (2) there are three consecutive iterations with- 
out a change in the design variables, or (3) a feasible 
direction is not found after 10 iterations. 

Design of Cylindrical Heat Pipes Embedded in 
a Honeycomb Core 

The cylindrical heat pipes embedded in the 
honeycomb core (fig. 6) were minimized with respect 
to mass, and the corrugations or channels were sized 
using the previously described procedures, subject to 
three heat transport constraints (sonic, entrainment, 
and axial wicking limits) (see appendix A, eqs. (Al) 
to (A3)). The design variables are the inner and 
outer radii of the channels (r, and r 2 , respectively) 
and the width of the grooves W. The heat pipes 
were first optimized assuming that the coolant man- 
ifolds were located at either end of the panel (20 ft 
apart). The panels were then optimized for a vari- 
able number of equally spaced manifolds N m , which 
accounted for the additional mass associated with 
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each manifold (0.065 lbm/ft 2 (ref. 7)). The working 
fluid chosen is toluene, and its properties at 300 r 
were obtained from reference 12. Hence, the^three 
design variables for this problem are: X\ - n, 
Xo = r 2 , and X 3 = W. The heat-pipe, heat- 
transport constraints (eqs. (Al) to (A3)) are related 
to the design variables through the cross-sectional ar- 
eas, A v and At, and through the expression for wick 
permeability. The resulting expressions for objective 
function and constraints are as follows. 

(p,A, + PsAs + *...A.»KNo. of heat pipea)L + M 
F[X) = bL 


gj(X) = (ri/ r 2) ~ 1-0 (geometric constraints) 

g 2 (X) = 2(r 2 +0.02) - 1.0 
g 3 (X) = YS(Q,n/Qaw) - 10 (wicking limit) 
g 4 (X) = FS (Qin/Qae) ~ 10 (entrainment limit) 
g 5 (X) = FS(Q/jv/Qas) - 10 (sonic limit) 

Upper and lower dimensional bounds (in inches) were 
also imposed as follows: 


0.05 < r\ < 10 


0.05 < r 2 < 1.0 


0.005 <W < 0.25 


No constraint was imposed to ensure an integer 
number of grooves; hence, one of the grooves could be 
nonstandard in shape. The limits Qaw, Qas, and 
Qae were calculated using equations (Al) to (A3) 
with 2 

A v = n r 2 
A t = rr {r\ - r?)/ 2 

r ,9 2 1 ( r 2 _ r l) 

A s = 7 T [(r 2 + 0.02 ) 2 - rl\ + 


for the axial-grooved design and 
At = rr(r 2 - r\) 

A fl = 7 r[(r 2 + 0.02) 2 -ri] 

for the corrugated composite screen wick. The per- 
meability K w is calculated for both wick designs as 


Ky 


2p2H -- f-12.956(D/VF) 3 

(. D + W ) 2 1 

+ 33.605(D/VF) 2 - 30.249 [D/W) + 24.0] 


where this expression represents a third-order curve 
fit of data presented in reference 13 and where 

D = r 2 — r\. . 

The heat input to a single heat pipe and the 

effective heat-pipe length are 


QlN ~ (9abs W 'eff) No of heat pi pe s 


and _ 

L e ff = Lj [2(N m + 1)] 

where L is the length of the panel (20 ft), <7 a bs 13 
the net heat absorbed by the panel per unit area 
(0.8 Btu/ft 2 -sec), and N m is the number of addi- 
tional manifolds. 

Design of Heat-Pipe Web-Core Sandwich 
Panel 

The heat-pipe web-core sandwich panel design 
(fig. 7) has a total of seven design variables. The 
design variables are the thickness of the face sheets 
t f , the height of the core h c , the web spacing S, 
the thickness of the core web t c , the thickness of the 
composite wick along the web walls t w , the depth 
of the channels D, and the width of the channels 
W. Hence, for this problem, the design variables are 

as follows: Xi = tt, X 2 = t c , X 3 = h c , X 4 = S, 
X 5 = t w , Xq = D, and X 7 = W. The constraints in- 
clude axial (wicking, entrainment, and sonic limits) 
and transverse (wicking limit through the depth o 
the panel) heat-pipe constraints and strength, buck- 
ling, fatigue, and minimum-gage constraints. The 
objective function to be minimized is the panel mass 
per unit planform area. The effect of additional 
coolant manifolds is also investigated. The steady- 
state, heat-pipe-limit equations are similar to equa- 
tions (Al) to (A3) of appendix A, with appropriate 
changes for At, A v , etc. A transverse wicking limit 
(through the depth of the panel) was also included 
in the analysis of the web-core design to ensure ade- 
quate liquid flow from the cooled to the heated face 
of the panel in the evaporator region (i.e., between 
the coolant manifolds). This amounts to a total 
of four heat-pipe constraints (three axial and one 
trsinsvcrsc] . 

In the structural analysis, it is assumed that the 
panel behaves as an orthotropic panel with panel 
stiffnesses calculated as functions of the design vari- 
ables using references 14 and 15. Details of the 
strength, buckling, and fatigue analyses are pre- 
sented in appendix B. The expression for the objec- 
tive function is given by the following equation: 

<o.A. + o,A a + P W A W ) (No. of heat pipes) L + M 

nx) = bi 
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The expressions for the normalized constraints are 
follows: 

31 (X) - FS(Q ain /Q aw ) - 

- 1.0 (axial wicking limit) 

92 W = FS [Q A in /Q AE ) - 

1.0 (axial entrainment 
limit) 

93 W = FS {Q ain /Q as ) - 

1.0 (axial sonic limit) 

94 W - FS{Qt I n / Qtw ) — 

1.0 (transverse wicking 
limit) 

9sW = FS(lV/£)) - 1.0 

(groove geometry 
constraint) 

9eW = FS((t c + 2 t w )/S) - 

1.0 (wick thickness 
constraint) 

97{X) = FS K,//ff ult ) - 1.0 

(ultimate strength 
of faces) 

28(^0 FS(r 2Xjmax /r u i t ) — 10 (shear strength of 

core) 

99p0 - FS(<7t c /t7 u ] t ) - 1.0 

(tensile strength 
of core) 

9ioW=FS(P I /P xo( ,)-i.o 

(overall panel buckling) 

9u (-Xj FS((<r Iimax 4 - P x bf 

(local buckling 

to 

J 

1 

t— 1 

o 

of faces) 

212 (X) = FS(P x b/A 3 )/<T x , lb 

(local buckling of 

- 1.0 

webs) 

2l3(^) = (FS(cr x> p)y/<7 u]t ) 

(burst failure of faces 

- 1.0 

in x-direction) 

214P0 — (FSfcry^p)^/^^) 

(burst failure of faces 

- 1.0 

in y-direction) 

2is(^) — (FS(t7 x ,p)c/<7 u ] t ) 

(burst failure of core 

- 1.0 

in x-direction) 

216(^) = (FS(c7y,p) c /cr ult ) 

(burst failure of core 

- 1.0 

in y-direction) 

2l7(^") = (o"i,max 

(fatigue failure of 

+ F > x /As)<J x ^A r T 
- 1.0 

faces) 

where 

Qain = Qin 
and 


QtIN = <?abs(S , £ eff ) 

Upper and lower dimensional bounds (in inches) were 
also imposed as follows: 

0.01 < t f < 1.0 
0.01 <t c < 1.0 
0.125 <h c < 10.0 


0.25 < 5 < 2.0 

0.005 <t w < 0.25 

0.005 < D < 0.25 

0.005 < W < 0.25 

Results and Discussion 

Cylindrical Heat Pipes Embedded in a 
Honeycomb Core 

Results for the optimization of cylindrical heat 
pipes embedded in a honeycomb core are summarized 
m table 1. The optimization problem can be eas- 
ily solved directly by noting that g 3 is the only con- 
straint that is a function of W via the expression for 
Ky. The optimum dimensions of the groove which 
minimize fluid flow resistance occur when W = r 2 - ri 
(ref. 13); hence, the problem can be reduced to two 
dimensions or design variables if the wicking limit is 
a governing condition of the design, which it happens 
to be for this problem. This is a valid assumption, 
since the optimum design occurs at the intersection 
of several constraint boundaries, which for this prob- 
lem are the entrainment-limit and the wicking-limit 
boundaries. Contours of constant values of the objec- 
tive function F (which for this problem is only the 
mass of the heat pipes) and constraint boundaries 
are plotted as a function of the design variables n 
and r 2 in figure 8. The feasible design space, with 
all constraints satisfied (g t < 0), is marked by the 
s a ed region. It is bounded by (a geometric con- 
straint on the outer diameter of the heat pipe) q 3 
(the wicking limit), and g A (the entrainment limit) 
ihe sonic limit is not a limit on the feasible design 
space. (If plotted, the sonic-limit boundary would 
be a vertical line located at r x = 0.07 in.) Contours 
ot constant objective function values are shown with 
the direction of increasing magnitude indicated The 
optimum design point indicated in figure 8 occurs 
at the intersection of the entrainment- and wicking- 
limit curves. This closed-form solution corresponds 

ri =°? 75 in > r \= 0343 = 0.068 in., and 

* 0 0145 Ibni/in . The problem was also solved 

using formal optimization methods and allowing W 
to be a free variable. The results after nine design 
iteration cycles from an initial infeasible design of 
n - 0.1 in., r 2 = 0.2 in., and W = 0.025 in. are 
n — 0.277 in., r 2 = 0.359 in., W = 0.0711 in. 
and F = 0.0149 lbm/in 2 . The total HPCP mass 
was 0.033 lbm/in after the face sheets were resized 
for structural constraints. The total panel mass 
(ret. 7) includes the heat-pipe mass and two mani- 
folds (0.065 lbm/ft each) but does not include the 
radiation system (1.49 lbm/ft 2 ), panel coolant (0.12 
lbm/ftZ), D-tube (0.16 lbm/ft 2 ), and difference in re- 
sized face-sheet masses. This mass compares with an 
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ACP mass of 0.02104 lbm/in 2 from reference 7 (to- 
tal RACP panel mass not considering the radiation 
system). The iteration history of the formal opti- 
mization procedure is shown in figure 8 by the ar- 
rows. The starting design is infeasible (violates con- 
straints 3 and 4); however, the process moves into 
the feasible region after three iterations and satisfies 
the first termination criterion after nine iterations 
(less than 0.1-percent change in the objective func- 
tion). Results from the direct solution and the formal 
optimization solution agree well with slight discrep- 
ancies caused by the approximate nature of the ex- 
pression for K w and by the fact that W f r 2 - n- 
If an additional constraint were included in the lor- 
mal optimization such that W = r 2 - r l> the re- 
sults would agree exactly with the direct solution. 
The simplicity of this problem enables a graphical 
solution for comparison; for practical design stud- 
ies where the number of design variables and con- 
straints is large, an exact solution is not possible^ 
Plots of the objective-function, design-variable, and 
constraint iteration histories are shown in figures 9 
to 11. As shown in figures 9 and 11, the initial design 
is infeasible with a low mass. The mass increases as 
the design moves toward the feasible design space 
and reaches a maximum after two cycles (fig. )• 
Figure 10 shows the iteration history of the design 
variables; the design variables converge after eight 
iterations. Designs are in the feasible design space 
after three iterations (fig. 11) and remain feasible 
with minor excursions, throughout the remainder of 
the optimization history. The desipi process pro- 
ceeds in the feasible region and minimizes the mass 
as previously shown in figure 8. This procedure was 
repeated for increasing numbers of additional man- 
ifolds and for the corrugated-screen composite wick 
design, and the results are presented in figure • 
A minimum-mass, axial-grooved, heat-pipe pane 
(0 0092 lbm/in 2 ) results for five equally spaced addi- 
tional coolant manifolds (3Vs ft apart); this amounts 
to a total panel mass of 0.0284 lbm/in 2 , compared 
with 0.033 lbm/in 2 if there are no additional man- 
ifolds. Results for the corrugated-screen design in- 
dicate a total HPCP panel mass of 0.0283 lbm/in 
with no additional manifolds and a minimum mass 
of 0.0276 lbm/in 2 with two additional manifolds. 

Heat-Pipe Web-Core Sandwich Panel 

Results for the screen-covered, grooved-wick de- 
sign (fig. 7(a)), starting from an initial design listed 
in table 2, indicate that with no additional mani- 
folds a minimum cooled-panel mass of 0.0279 lbm/in 

reSU The total HPCP mass is the sum of the total ACP 
mass of reference 7 (4.52 lbm/ft 2 ) and the heat-pipe 


web-core panel mass with any additional coolant 
manifolds (0.065 lbm/ft 2 each) minus the face sheets 
(1.20 lbm/ft 2 ), D-tubes (0.16 lbm/ft ), honeycomb 
(0.29 lbm/ft 2 ), radiation system (149 lbm/ft ), 
panel coolant mass (0.12 lbm/ft 2 ), and adhesive 
(0 4 lbm/ft 2 ). A minimum mass of 0.0221 lbm/in 
is achieved by adding three additional manifolds^ 
Results for the corrugated-screen composite wick 
also listed in table 2 indicate that with no addi- 
tional manifolds, a minimum structural panel mass 
of 0.0161 lbm/in 2 results, which corresponds to a to- 
tal HPCP panel mass of 0.0221 lbm/in . The op- 
timum design with respect to additional manifolds 
is 0 0157 lbm/in 2 with one additional manifold and 
corresponds to a total HPCP panel mass of 0.0216 
lbm/in 2 . Again, this is compared with an ACP panel 
mass (ref. 7) of 0.02104 lbm/in 2 (total RACP panel 
mass excluding the radiation system). 

Iteration histories for a typical optimization are 
shown in figures 13 to 15 for a corrugated-composite 
screen wick design with no additional manifolds. e 
variable histories (fig. 14) were all normalized to 
their initial dimensions. Convergence to a local min- 
imum occurs after 18 iterations. The initial start- 
ing design is infeasible (violates axial-wicking and o- 
cal web-buckling constraints gi and g 12 , respectively, 
fig. 15). A feasible design is reached after two it- 
erations (F = 0.0174 lbm/in 2 ), and panel mass is 
reduced to a minimum value of F - 0.016 lbm/in 
after 18 iterations (fig. 13). 

The individual CPU times needed to achieve the 
minimum HPCP designs using formal optimization 
methods are about 1 minute with a CDC Cyber 173 
computer. Hence, if new designs are required because 
of changes in loading conditions, constraints, etc., op- 
timum configurations can be determined quickly. Al- 
though formal optimization techniques may appear 
to be useful in a black-box mode, it is necessary 
for someone familiar with each aspect of the engi- 
neering design to monitor results of the optimization 
process. Such engineering experience is necessary to 
make critical decisions when problems arise m the 
optimization because of 

1. a failure to satisfy specific constraints and pro- 
duce a feasible design 

2. a failure to include a relevant design constraint 

3. inadequate analysis methods based on changes 
which occur during the design 

4. any number of other scenarios which might pro- 
duce a design which does not appear to be 
realistic 


8 


Concluding Remarks 

Results of this study indicate that formal op- 
timization methods can be used to design heat- 
pipe systems to perform several functions efficiently. 
Several structurally efficient heat-pipe panels were 
designed and optimized with respect to mass and 
subjected to heat-transport, strength, buckling, fa- 
tigue, and mmimum-gage constraints. Optimum 
panel designs were also obtained for various num- 
bers of additional coolant manifolds where the ther- 
mal but not structural effect of additional manifolds 
was considered. A nonmtegrated design— cylindrical 
heat pipes embedded in a honeycomb sandwich- 
with three design variables and five constraints, was 
optimized using an optimizer, CONMIN, for con- 
strained function minimization. Optimum panel 
masses for a screen-covered, axial-grooved wick de- 

3 00284 ' f m / in T th n ° additio -l manifolds 
and 0.0284 lbm/m 2 with five additional manifolds. 

the corrugated-screen wick design, the optimum 
panel masses were 0.0283 and 0.0276 lbm/in 2 with 

foLk ddltl ° na ! m . anifo,ds and additional mani- 
folds respectively. The accuracy of the optimum 

fo m ll! ng LONMIN was checked w >th a closed- 
orm solution for the cylindrical heat pipe, and re- 
sults agree well. Results for an integrated design— 


heat-pipe web-core sandwich panel-were consis- 
tently lighter than for the discrete tube design (heat 
pipes embedded in honeycomb sandwich) by ap- 
proxmiatdy 20 percent. Results for screen-covered^ 

0 n^i ed iK aC< ^ S 2 eetS ,_ W,ck design were 0-0279 and 
W< th "° addit i° n al manifolds and 
ee additional manifolds, respectively. Results for a 

wick design were 0.022 

“ d °,°, 2 ! 6 lbl f/ ln with no additional manifolds and 
one additional manifold, respectively. The web-core 
panel design had a total of 7 design variables and 
17 constraints; hence, formal optimization methods 
were necessary to find a minimum-mass design Fea- 
sible designs were found after several design itera- 

abo u ; 2nV dGSignS Were rea * hed ^ 

about 2° iterations and required about 1 minute of 
CPU time on a CDC® Cyber 173 computer. 

The study also verified the feasibility of using 
heat pipes in conjunction with active cooling to cool 
the structure of hypersonic cruise aircraft. The 
resulting designs offer a mass-competitive alternative 
a redundant, completely actively cooled design. 


NASA Langley Research Center 
Hampton, Virginia 23665-5225 
August 14, 1987 
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Appendix A 

Heat-Pipe Operational Limits 

Heat-pipe operation is analyzed assuming steady- 
stJe, one-dimensional, incompress, ble, — « 
flow with complete pressure recovery m the s vapo 
region Liquid flow in the wick is described by 
Darcy’s law and a constant heat flux is assumed 
over Y the panel surface. Three fluid-dynamic heat 
tr^sfer limits, the wicking limit (capillary pumping 
limit), the entrainment limit, and the sonic limit, 
were considered in the present analysis. 

Wicking Limit 

The wicking limit is the heat-transport limit that 
occurs when the capillary pumping action of t 
is iust sufficient to balance the gravitational pressur 
head of the liquid and the inertial and frictional losses 
of the fluid. The expressions for the wicking lim 
obtained as follows using the method of reference . 

Q aw = 2 a t ! {r p [{livL e ff/gcA t K w \pt) 

+ (8 ntL e ff/gc^ r tPv^)\} 

where at is the liquid surface tension, p v and ^ 
are the vapor and liquid mass densities r p is the 
^ore radius of the wick, A, is the liquid-flow cross- 
sectional area, A is the latent heat of vaporization of 
the working fluid, r v is the 

u , nnr o Dace L„ff is the effective length of the heat 
pipe u v and M are the vapor and liquid viscosities, 
AT P is the permeability of the wick, and g c is the 
pavLt,aTcon,ta n , Since coolant manifolds are 


located at either end of the panel, it is assumed 
that the effect of vehicle acceleration forces can be 

neglected. 

Entrainment Limit 

The entrainment limit occurs when the vapor 
velocity is large enough so that the drag force exer 
on the returning liquid condensate ^erco^the 
surface-tension forces of the liquid and liquid droplets 
are entrained in the vapor space and swep .away i 
the opposite direction. This causes liquid depletion 
in the heated area and the eventual drying out of the 
heat pipe. The entrainment limit is assumed to occur 
when the Weber number (ratio of mertial vapor force 
to liquid surface-tension force) is equal to unity. T 
entrainment limit is expressed as 

Qae = A v (gcpv°i>? / 2r p) / 


Sonic Limit 

The sonic limit occurs when the vapor flow veloc- 
ity reaches a sonic condition. The heat transported 
is thus limited by this choked flow condition and is 
expressed as 


Qas = 


AvPv Ado 

V / 2(l + 'l) 


(A3) 


where A v is the cross-sectional area of the vapor 
space, p v is the vapor pressure, l is the ratio ot 
specific heats, and a G is the sonic velocity of the 

vapor. 
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Appendix B 

Structural Analysis of Heat-Pipe Web-Core 
Sandwich Panel 

In the structural analysis the panel is assumed to 
behave as an orthotropic panel with stiffness prop- 
erties calculated as functions of the design variables 
using references 14 and 15. Equations derived in ref- 
erence 17 are used to calculate maximum bending 
moment and shear load for a simply supported panel 
(2 ft x 2 ft, dimensions of panel between rings and 
stiffeners) subject to in-plane and uniform transverse 
loads as shown below: 


M x .max — qb 2 02 (Bl) 

Qx,max — qb(3§ (B2) 

where and /% are functions of all the design 
variables shown in figure 7 except t w (ref. 17), where 
q is the applied uniform transverse load, and where 
b is the width of the panel. The panel bending stress 
is then calculated using simple beam theory as 


0a;, max — ^x,ma xy/I x (B3) 

where y is the distance from the neutral axis of the 
panel to the midplane of the faces (y = (tf + h c /2)) 
and where I x is the moment of inertia of the panel 
about the neutral axis per unit panel width. The 
maximum shear stress ^zx^msx is calculated from 
reference 17 to occur at the center of the core 
and is 

r zx, max = V Q x ,max/ (2 I x t c ) (B4) 

where V is the static moment of the cross-sectional 
area above the neutral axis. 

The effect of internal pressure (p v = 39.9 psi for 
toluene at 300°F) was included using reference 18 
and assuming that the face sheets and core web 
(fig. 7) act as simply supported plates. The expres- 
sions for stresses in the plates are 


°x,p = |e P „/ 7T 2 t 2 (i/t 2 + l/w 2 ) 2 | (i/e 2 + v /w 2 ) (B5) 


7 y,p 


= {6p„/ pt 2 (i/e 2 + 1/™ 2 ) 2 ] } (v/e 2 + i/w 2 ) (B6) 


where the local width w (used in eqs. (B5) and (B6)) 
of the face-sheet plate is S, and the local width of the 
web plate is h c . The local length f is the length of the 
panel between rings (2 ft). The internal pressure is 
Pv, t is the plate thickness, l and w are the length and 


width of the plate, respectively, and v is Poisson’s 
ratio. Hence, the axial stress in the faces of the plate 
is 

°xj = ^i.max + °x,p + Pxb/ (B7) 

where P x is the applied in-plane load per unit width, 
-4 a is the total cross-sectional area of the panel which 
sustains load, and b is the panel width. The faces and 
walls were also sized to withstand, without failure, 
a burst pressure of four times the actual internal 
pressure. The tensile stresses in the webs due to 
internal pressure are 

otc = Pv(S/t c ) (B8) 

Overall panel buckling load is calculated using the 
following equation from reference 17: 

P x,ob = ($o*- 2 Ar) / [b 2 (m 2 fc 2 /^ 2 )] (B9) 

where m is the number of half waves in the longitu- 
dinal direction and is varied to produce a minimum 
value of P x ob for buckling, and where <f> 0 is a buck- 
ling parameter defined in reference 17. 

Local buckling of the face sheets and webs was 
calculated using reference 19, assuming simply sup- 
ported plates, and results in 

a xjlb = Z-2V(t/w) 2 [E/(\ - v 2 )} (BIO) 

where t is the thickness of the plate and tv is the plate 
width ( S for face sheet and h c for webs). 

The panels were designed to sustain design limit 
loads without yielding and ultimate loads without 
failure. A modified Wheeler model was used for 
fatigue analysis, similar to that of reference 20, and 
material property values for aluminum alloy 
AL 2219-T81 were obtained from references 7 and 8. 
For values of the maximum stress ratio R = — 1 
(i? = /min //max), a, stress- intensity factor 

K c = 105 ksi- Vnu, and an initial flaw size of 0.005 in., 
the following equation was used to predict crack 
growth: 


da/dN = 3.8637 x 10 -6 (AA:) 2 2044 / 

[{l-R)K e -AK) (Bll) 

where 

AK — 1.12 (/max — /min) ^ (B12) 

where a is the crack length, Q is the flaw shape 
parameter (ref. 20), and / max and / min are the 
maximum and minimum cyclic stresses, respectively. 
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Plots of cycles to failure versus stress for various 
face-sheet thicknesses were used to determine an 
allowable stress as a function of face-sheet thickness 
(based on 20000 cycles to failure) to be used as the 
fatigue constraint. Interpolated values using a third- 
order curve fit results in the expression 


°x,FAT = 0-4608 x 10 8 («/) 3 - 0.8906 

X 10 7 (t/) 2 + 0.6363 x 10 6 (f/) 

+ 0.082727 x 10 3 (B13) 

for an allowable stress due to fatigue. 
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Table 1. Results of Cylindrical Heat Pipes Embedded in a Honeycomb Sandw.ch 


Cylindrical heat pipes 
in 

honeycomb 


Screen- covered grooves 


Corrugated-composite 

screen wick 


Initial 

design, 

in. 

Final 

design, 

in. 

^ Final cooled- 
panel mass, 
lbm/in 2 

ri = 0.2 

f‘2 ~ 0.3 
W = 0.05 

r 1 - 0.277 
r 2 = 0.351 
W = 0.0711 

0.033 

r\ — 0.2 
r2 = 0.3 
W = 0.05 1 

r, = 0.1665 
r 2 = 0.2187 
W - 0.0518 

0.0284 

r\ = 0.2 
T 2 0.3 
W = 0.05 

f j = 0.2554 
r 2 = 0.3222 
W - 0.0668 

0.0283 

n = 0.2 
r 2 = 0.3 
W = 0.05 

I n = 0. 1722 

r 2 = 0.2318 
W = 0.0593 

0.0276 


Radiation system not included. 


"V 


Design variables 

rj, r 2 , W 


Screen 


• 3 -h M -02 



Table 2. Results of Web-Core Sandwich Panel 
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(b) Parametric and trade studies. 

Figure 2. Design methodology used for radiantly and actively cooled panel (RACP (ref. 7)) 
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(a) bcreen-covered axial grooves. 


--J-WjWj. 



Working fluid 
passages 


(b) Corrugated sandwich screen wick. 

Design variables: t f , t c , h c> s, t w> D, W 


Figure 7. Heat-pipe web-core sandwich panel designs 
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Feasible design space 
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’ r = 0.0009055 in.; working fluid is toluene. 
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Figure 10. Design- variable iteration 


history (axial-grooved heat pipe with no additional 


coolant manifolds). 






1.6 


Normalization 

Variable factor, in. 



0.0300 

.0200 

1.0000 

.5000 


Figure 14. Design-variable iteration histones (web-core 
additional coolant manifolds). 


sandwich with corrugated-screen composite wick with 


no 


Normalization 



Figure 14. Concluded. 
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Figure 15. Continued. 
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Figure 15. Continued. 
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Figure 15. Concluded 
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